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Abstract 

We review some of the recent results on two-interacting particles (TIP) in low-dimensional disordered quantum 
models. Special attention is given to the mapping of the problem onto random band matrices. In particular, 
we construct two simple, seemingly closely related examples for which an analogous mapping leads to incorrect 
results. We briefly discuss possible reasons for this discrepancy based on the physical differences between the TIP 
problem and our examples. 

PACS: 71.55.Jv, 72.15.Rn, 71.30.-Hh 

Keywords: disordered systems, electron-electron interactions, TIP. 



1. Introduction 

Until 1994, the theoretical and experimental re- 
search on transport in disordered systems clearly 
supported the scaling hypothesis of localization 
for non- interacting electrons [y||]. The systems 
studied usually fell into the predicted universal- 
ity classes, and, if they didn't, then they could 
be shown not to be generic |^. However, real 
electrons of course interact H) , and their interac- 
tion is of relevance for the transport properties 
of disordered systems |^,|^ , especially in 2D and 
ID where screening H is less efficient than in 3D. 
The influence of weak interactions has been in- 
vestigated extensively using perturbation theory 
and the perturbative renormalization group (RG) 
ItIB]. One of the key results is that the lower 
critical dimension of the MIT is d^ = 2 as it is 
for non-interacting electrons. The application of 
the perturbative RG in ID p-[l3|] has lead to the 
prediction that all thermodynamic states remain 
localized in the presence of repulsive many-body 
interactions. 

Due to the persistent current problem [ |l4|-|l9| 
and new experiments on 2D electron systems |20|- 



E2J which show striking signatures of a metal- 
insulator transition, these theoretical considera- 
tions received a lot of renewed attention. In order 
to theoretically study the effects of the interplay 
between disorder and interactions, one should in 
principle solve a problem with an exponentially 
growing number of states in the Hilbert space 
with increasing system size. At present, this can 
be achieved only for a few particles in ID ||2^ 
P6[ and very few particles in 2D [27|~|3l|]. How- 
ever, in 1994 Shepelyansky ||3^ , |33 | proposed to 
simply look at two interacting particles (TIP) in 
a random environment. In particular, he sug- 
gested that the two particles would form pairs 
even for repulsive interactions such that the TIP 
pairs would have a larger localization length than 
the two single particles (SP) separately. Thus 
the interaction would lead to an enhanced possi- 
bility of transport p4| . The perhaps even more 
surprising part of the prediction is that the TIP 
pairs will have a localization length A2 such that 
at pair energy E — Q 



A2 ex f/^Ai' 



(1) 
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where U represents the onsite interaction strength 
and Ai is the SP localization length. Since Ai oc 
105/VK2 H-Hl in ID, this implies large values 



of A2 for small disorders W. 

The first numerical studies devoted to the TIP 
problem used the transfer-matrix method (TMM) 
to investigate the proposed enhancement of 
the pair localization length A2 [||,|§. The TMM 
of [Q contained an additional artificial infinitely- 
long-ranged interaction that tends to mask the 
onsite interaction Q. The TMM of |40) avoids 
this problem, but is restricted to small system 
sizes and results for localization lengths of A2 ~ 
300 had been deduced on systems of size M = 
100. Therefore, two of us studied the TIP prob- 
lem by a different TMM |E3] at large system size 
M > 300 and found that (i) the enhancement 
A2/A1 decreases with increasing M, (ii) the be- 
havior of A2 for L'^ = is equal to Ai in the limit 
M — > cx) only, and (iii) for [/ 7^ the enhance- 
ment A2/A1 also vanishes completely in this limit. 
Consequently, we concluded g] that the TMM 
applied to the TIP problem in ID measures an 
enhancement of the localisation length which is 
due to the finiteness of the systems considered. 
Although Ref. [|| has been criticized [||,|3, 
we emphasize that subsequent publications have 
shown I^SJ-^ that there are no variants of TMM 
that reproduce Eq. ( [l|). Furthermore, in a later 
numerical approach [[48| , based on Green func- 
tion methods. Song and v. Oppen argue that our 
extrapolations for M ^ 00 were off by « 11% 
only, whereas the original TMM of |Q deviated 
by about a factor of 3 |4^ . Thus while our criti- 
cized TMM results are valid, various other numer- 
ical investigations by other groups [Eq-M as well 
as ourselves |58~3l| convincingly demonstrated 
some enhancement. The reason for the failure of 
the TMM approach of p^Ji^ has been explained 
by Song and v. Oppen ||4^ by arguing that the 
TMM measures a localization length A/ < A2 due 
to the cigar-shape geometry p2] of the TIP states. 

Reliable numerical approaches to the TIP prob- 
lem are nowadays based on the computation of 
the decay of the Green function |8|J5|j5^,||,|| . 
Other direct numerical approaches to the TIP 
problem have been based on the time evolution of 
wave packets ||3^,^,^,^, exact diagonalization 
164 , variants of level statistics] 4^ 51[ and analy- 
sis of multifractal properties [£4,^|, perturbative 
methods [|64|j65[| and mappings to effective models 



-|70[. In these investigations an enhancement 
of A2 compared to Ai has been found as remarked 
above but the quantitative results tend to differ 
both from the analytical prediction in Eq. (|l|), 
and, albeit less, from each other. Furthermore, a 
check of the functional dependence of A2 on Ai is 
numerically very expensive since it requires very 
large system sizes M ::^ A2 ^ Ai. Extensions 
of the original arguments have been proposed for 
TIP in 2D l2^j3^,|6l|, |69|,[7l| ,[72| and 3D |7|, for TIP 
close to a Fermi sea |74], and for long-range in- 
teractions in ID [|l]j4|]5|. 

The basic idea leading to the prediction ffl) 
is based on looking at the interaction matrix 
element between two eigenstates tpki = 4'kipi 
and ^pnm = ipn^'m of the non- interacting sys- 
Here ipk,iph^pn,^pm denote SP eigen- 



32 3 



tem 

states localized with localization length Ai around 
sites k,l,n,m. For an onsite interaction iTq] 
V^.^njiUj^ (with rijT denoting the number 
operator at site j for spin r) only states with 
\k-l\ < Ai, \n-m\ < Ai, \k-n\ < Ai, \l-m\ < Ai 
will give significant contributions to the interac- 
tion matrix element 



(^pkl\U\^p„ 



N 



c^Ev'I(j')W(^>"(j>™(j')- 



(2) 



i=i 



These conditions are illustrated in Fig. |l]. If one 
assumes l33p3] that the SP state is given as 



i^kij) oc 



/AI 



exp 



\3 - k\ 
Ai 



^OU) 



(3) 



with 6{j) a random phase, one finds [^ that the 
typical interaction matrix element has a magni- 
tude of 



u ex A 



-3/2 



(4) 



since it is the sum of Ai random contributions 
of magnitude A^ . Shepelyansky next calculated 
the decay rate F of a non-interacting eigenstate 
by means of Fermi's golden rule F ~ U^/Xit 
P^ , p3| , |64| . Since the typical hopping distance 
is of the order of Ai the diffusion constant is 
D ^ U^Xi/t. Within a time t the particle pair 

states. Diffusion stops 



visits N - UXl^^t 



-1/2^1/2 



Tnm Tn Tm 




«c > 

Figure 1. Schematic picture of the TIP argu- 
ments of Ref. [^. The two-particle state ipki 
(left solid and dashed exponentials indicate the 
envelopes of the constituents ipk and ipl) is local- 
ized within a distance Ai from the two-particle 
state ipnm (right solid and dashed curves). The 
resulting overlap-matrix element u — (U) = 
(tpkilUltpnm) leads to a longer decay length A2 
for the TIP state as explained in the text. This 
effect can be visualized as an effective reduction 
(thick short-dashed line) of the original disorder 
potential (thin short-dashed line). 



when the level spacing of the visited states is of 
the order of the frequency resolution 1/r. This 
determines the cut-off time r* and the corre- 
sponding pair-localization length is obtained as 
A2 ~ \JDt* ~ {U/t)'^\\ in agreement with Eq. 
(El). Applicability of Fermi's golden rule requires 
r » t/\\ which is equivalent to U'^Xi/t^ > 1. 
This is exactly the condition for an enhance- 
ment of A2 compared to Ai. Alternatively, the 
model may be mapped to a random-matrix model 
(RMM) with entries chosen according to Eq. (|j) 



making the arithmetic average unsuitable to char- 
acterize the typical value. In Ref. [Q we have 
paid special attention to the exact dependence of 
u on Ai and system size. To this end, we di- 
agonalized the ID Anderson model for a given 
M and W and computed u by averaging over 
all suitable states and many disorder configura- 
tions. We showed that due to the strongly non- 
Gaussian distribution of u, one should rather use 
the logarithmic average than the arithmetic av- 
erage as the typical value for the computation 
of m(Ai). But whereas the arithmetic average 
||54|] gives Uabs oc Aj"^'^, the typical value obeys 



Utyp OC 



A: 



-1.95 



Following the arguments above. 



this would imply A2 oc A}"'^, i.e., a very small en- 
hancement. We emphasize that this result does 
not mean that there is no enhancement of the 
localization length. Rather, the results of Ref. 
fZ^ ] indicate that the arguments of Ref. |Q cap- 
ture the physics, but only in a somewhat sim- 
plified form. One step towards a better agree- 
ment between the analytical and the numerical 
approaches is to take into account the energy de- 
nominators in the computation of u, e.g., to con- 
sider only interaction matrix elements for states 
whose energy spacings are of the order of U or 
smaller ^^. In this case we find that there is 
a slight decrease in the value of the typical ex- 
ponent and correspondingly a slight increase in 
TIP delocalization yielding A2 oc A}-^^°^ This 
suggests that higher orders in perturbation the- 
ory than the first order RMM approach ||3|] are 
important. Furthermore, the exponent 1.4 ±0.2 
is in reasonable agreement with previous results 
in the literature P,P,p|-p5|,p7|J5S|J6^ p|,p9|,ra . 



2. Numerical results for the random- 
matrix model of TIP 

The arguments presented in the last section 
are of a qualitatively nature and Eq. (|l|) must 
be checked for quantitative accuracy. Even be- 
fore testing (^, it is already worthwhile to check 
the validity of (P and the subsequent arguments 
or the RMM approach ||,^,|||. In Ref. [|o), it 
had been shown that the assumption of a Gaus- 
sian distribution of the matrix elements u was 
oversimplified. The distribution showed long tails 



3. RMM approach for toy models 

In this section we show that a naive application 
the RMM approach may give qualitatively incor- 
rect results even if the RMM contains the correct 
dependence of the matrix elements on the disor- 
der strength. To this end we consider two toy 
models which seem to be closely related to the 
TIP problem. For these models, viz. Anderson 
models of localization with additional perturbing 
random potentials, we show that mapping onto 
RMMs and estimating the localization length by 



Fermi's golden rule leads to an incorrect enhance- 
ment of the localization length. 

3.1. 2D Anderson model with perturba- 
tion on a line 

The first example is set up to lead to the same 
RMM as the TIP problem. It consists of the usual 
2D Anderson model perturbed by an additional 
weak random potential of strength U at the diag- 
onal X = y in real space. Since the perturbation 
increases the width of the disorder distribution 
at the diagonal we expect it to decrease the lo- 
calization length. We map the model onto an 
RMM following the arguments for the TIP prob- 
lem sketched in Sec. |l|. Again, the eigenstates of 
the unperturbed system are localized with a lo- 
calization length Ai and approximately given by 



tp„{x,y) r^ — exp 



\r - r„ 
Ai 



iOn{r) 



(5) 



where r = {x,y) is the coordinate vector of the 
particle and 9 is again assumed to be a random 
phase. The Hamiltonian of the 2D perturbed An- 
derson model differs from the TIP Hamiltonian 
in two points: (i) the diagonal elements (given 
by the random potential) are independent ran- 
dom numbers instead of being correlated as in 
the TIP problem and (ii) the perturbing poten- 
tial U{x, x) £ [—U, U] at each diagonal site is ran- 
dom instead of having a definite sign and modu- 
lus U as in the TIP problem. However, none of 
these points enters the mapping procedure out- 
lined in Sec. n^. Thus, we find that the pertur- 
bation couples each state close to the diagonal 
{\xn — yn\ < Ai) to O(A^) other such states. 
The interaction matrix element is again a sum 
of 0(Ai) terms of magnitude U/\1 and random 
phases and as before u ^ UXi . Consequently, 
our toy model is mapped onto exactly the same 
RMM as TIP in a random potential. Therefore, 
the resulting localization length along the diago- 
nal is also given by Eq. (|]) . We thus arrive at the 
surprising conclusion, that adding a weak random 
potential at the diagonal of a 2D Anderson model 
leads to an enormous enhancement of the local- 
ization length along this diagonal, in contradic- 
tion to the expectation expressed above, viz. that 
increasing disorder leads to stronger localization. 
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Figure 2. Dependence of Ai(M) on disorder 
W for the 2D Anderson model at E = for 
M = 10, 25, 30, 35 and 50 indicated by increasing 
symbol size. We use the M = 50 data, empha- 
sized by the solid line, as finite-size estimate of 
Ai. 



As for the TIP case [^ we now numerically 

3/2 

check whether the relation u ~ UXi between 
the coupling matrix element u and the localiza- 
tion length Ai of the unperturbed system is cor- 
rectly described by the RMM. Since in 2D a sim- 
ple analytic formula for the dependence of Ai on 
the disorder W does not exist, we first compute 
estimates Ai(M) for quasi-lD strips of finite strip 
width M with 1% accuracy by TMM. In Fig. |, 
we show data of Ai (M) as a function oiW. In the 
following, we take Ai (50) to compute the coupling 
matrix elements. 

Next, we calculate both the arithmetic aver- 
age Mabs — {\u\) and the logarithmic average 
Mtyp — exp[(log(|u|))] for different values of W 
and various M x M squares. Disorder averaging 
is over 20 samples and we study Uabs and Mtyp as 
functions of Ai(M). We emphasize that instead 
of the well-known extrapolations of Ai(M) to in- 
finite system size by means of FSS 0, we take 
the finite-size approximants Ai(Af) on purpose, 
since we compute A2 also for comparable finite 
sizes only. 

As for the TIP model (76) the distribution 
Po{u) of the (off-diagonal) coupling matrix ele- 
ments is strongly non-Gaussian, suggesting that 
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Figure 3. Dependence of Uabs (squares) and utyp 
(circles) on Ai(M) for the perturbed 2D Anderson 
model with [/ = 1 and Af = 10, 25, 30 and 35 
indicated by increasing symbol size. The solid 
lines represent the power laws Wabs ~ A|j~ ' and 



Wtyp 



AT 



Wtyp rather than Uabs is the relevant quantity. 
The results for Uabs and wtyp are presented in 
Fig. H The dependence of Uabs on Ai(M) for 
2 < Ai(M) < 12 follows Wabs oc Ai(M)-i-6±o.i 
in agreement with the RMM value of 3/2 and 
with Ref. JT^]. Furthermore, here we also have 
wtyp oc Ai(M)~^-^^°'^. We note that the change 
of the slopes of Uabs and Utyp at Ai(M) ^ M/2 is 
entirely due to the finite sample sizes ^^ . 

Consequently, in contrast to the TIP problem 
the RMM model for the 2D perturbed Anderson 
model of localization contains the correct depen- 
dence of the coupling matrix elements on the lo- 
calization length of the unperturbed system, but 
still it leads to an incorrect enhancement of the 
localization length along the diagonal. 

3.2. ID Anderson model with perturba- 
tion 

An even more striking contradiction can be ob- 
tained for an ID Anderson model of localization. 
The eigenstates are again given by Eq. (y) with 
Ai known from second order perturbation the- 
ory ||35|-|37[| and numerical calculations IUJI^ to 
vary as Ai ~ t'^ /W"^ for small disorder. We now 
add a weak random potential of strength U at all 



sites. Since the result is obviously an ID Ander- 
son model with a slightly higher disorder strength 
Wu > W the localization length will be reduced, 

A(c/) ^ tywi 

The mapping onto an RMM can be performed 
in complete analogy to the TIP problem and the 
2D Anderson model discussed above. The per- 
turbing potential leads to transitions between the 
unperturbed eigenstates ipn- Each such state is 
now coupled to 0(Ai) other states by coupling 
matrix elements {ipnlUltpn') with magnitude u '^ 

— 1/2 

UXi since we sum over Ai contributions with 
magnitude U/Xi and supposedly random phases. 

The application of Fermi's golden rule in this 
ID case leads to a diffusion constant D ~ U^Xf/t. 
The number of states visited within a time r is 
N ^ f/Aii~^/^r^/^. Again, diffusion stops at a 
time T* when the level spacing of the states vis- 
ited equals the frequency resolution. This gives 
T* ~ U'^Xl/t^. The locahzation le ngth A of the 
perturbed system thus reads A ^ \J Dt* ~ f/^A^ 
as in Eq. (|l|), in contradiction to the correct re- 
sult. 

Again we numerically check the relation be- 
tween Uabs and utyp and the unperturbed local- 
ization length Ai. In Fig. ^, we show results 
obtained for chains with various lengths and 50 
disorder configurations for each W . Ai is com- 
puted by TMM. For 10 < Ai < 250, Uabs varies 
as A^""*^ °'^° as predicted above. Utyp varies as 
^-o.59±o.io gQ^j^ variations are compatible with 
the RMM value of 1/2 for the exponent. Again 
we need Ai < M jl in order to suppress finite size 
effects. 

Consequently, although the RMM model for 
the ID perturbed Anderson model of localization 
contains the correct dependence of the coupling 
matrix elements on the localization length of the 
unperturbed system, it still leads to an incorrect 
enhancement of the localization length. 

4. Application of the block-scaling picture 
to toy models 

Let us now discuss the relation of these results 
to Imry's block-scaling picture (BSP) p^ , |69[ for 
the TIP problem. In this approach one considers 
blocks of linear size Ai and calculates the dimen- 
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Figure 4. Dependence of Mabs (squares) and utyp 
(circles) on Ai for the perturbed ID Anderson 
model with f7 = 1 and M = 200, 300, 500 and 
800 indicated by increasing symbol size. The solid 
lines represent the power laws Uabs ^ ^i and 
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0.59 



sionless pair conductance on that scale, 



52 



u 

A2' 



(6) 



where u represents the typical interaction- 
induced coupling matrix element between states 
in neighboring blocks and A ^ t/\\ is the level 
spacing within the block. If the typical coupling 
matrix element depends on Ai as u ~ C/A^" the 
pair conductance obeys 



92 - [U/ifXt 



(7) 



For the 2D Anderson model with perturbation 
considered above, the BSP can be applied anal- 
ogously. Again, we consider blocks of linear size 
Ai and compute the typical perturbation-induced 
matrix elements between these blocks. We then 
find that according to the BSP the conductance of 
a 2D Anderson model with additional weak per- 
turbing potential along the diagonal is given by 
Eq. (^. Using a = 1.5 ± 0.1 as obtained above 
from the numerical data for Uabs and utyp, we 
then have 52 ~ {U/t)'^Xi. Thus we conclude that 
the BSP does not work for our 2D toy model, be- 
cause it yields the same unphysical result as the 
RMM approach of section 3.1, 



Let us also apply the BSP to the ID toy exam- 
ple. The level spacing in a ID block of size Ai is 
A ~ t/Xi, and the coupling matrix element be- 

— 1/2 

tween states in neighboring blocks is m ~ UXi . 
Thus, the conductance of the perturbed system 
on a scale Ai is obtained as (72 ~ (C//t)^Ai. For 
large Ai this again contradicts the correct result, 
viz. a decrease of the conductance compared to 
the unperturbed system. Thus, the BSP applied 
to the two toy models introduced in the present 
work gives the same qualitatively incorrect results 
for the localization properties as the RMM. This 
is not surprising since the only ingredients of the 
BSP are the intra-block level spacing A ^^ i/Af 
and the inter-block coupling matrix elements u 
which also enter the RMM. 



5. Conclusions 

We have presented two toy models which seem 
to be closely related to the TIP problem. For 
these toy models the usual analytical arguments 
given to support the delocalization of TIP, viz. 
the RMM and the BSP do not work. However, 
the large-scale numerical simulations p8| , p0|j5^ , 
py,0,^ have convincingly shown that an en- 
hancement of the two-particle localization length 
due to the interaction exists, even though the de- 
tailed results are more complicated than the orig- 
inal prediction (p. This leads, of course, to the 
question, under which conditions the RMM map- 
ping and the BSP give the correct result, at least 
qualitatively. While a general answer to this ques- 
tion is not known, it has been suggested |Q that 
the difference between the TIP and our toy mod- 
els is an additional symmetry in the TIP problem. 

In summary, the two examples suggest that ad- 
ditional physical insight is needed before apply- 
ing the RMM. In addition, we expect that taking 
into account the energies of the states as in Ref. 
pQ l for TIP will result in a reduced enhancement, 
i.e., a smaller value of a, in the analytical predic- 
tions. This will in turn give a better agreement 
with the numerically determined dependence of 
the TIP localization length on Ai. 
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